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^ ; Abstract 

We present a next-to-leading order (NLO) computation of the full set of polarized 
and unpolarized electroweak semi-inclusive DIS (SIDIS) structure functions, whose 
knowledge is crucial for a precise extraction of polarized parton distributions. We 
focus on the phenomenology of the polarized structure functions for the kinematical 
conditions that could be reached in an Electron-lon-Collider. 

We show that the NLO corrections are sizeable, particularly in the small-x range. 
We test the sensitivity of these structure functions on certain quark distributions and 
compare it to the situation of inclusive DIS and electromagnetic SlDlS. 
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1 Introduction 



Understanding how the nucleon spin is composed of the angular momenta and spins of its con- 
stituents (quarks and gluons) has been a defining question in hadron structure for a long time. 
Since it was found that little of the proton spin is carried by the quarks and anti-quarks spins, sev- 
eral experiments have measured with increasing precision observables which are sensitive to quark 
and gluon polarizations in the nucleon. Such experimental progress was matched by advancements 
in theoretical precision and phenomenological analyses of data [1]. 

The spin structure of a nucleon can be described by the (anti) quark and gluon polarized parton 
distribution functions (pPDFs), defined by 

Af,{x,Q') ^ f;{x,Q') - f-{x,Q'), (1) 

where fj~{x, Q^) {f[{x, Q"^)) denotes the distribution of a parton j with positive (negative) helicity 
in a nucleon with positive helicity, as a function of momentum fraction x and scale Q. Its first 
momentum, that is to say the integral Afj(Q'^) = Afj(x,Q'^)dx, directly measures the spin 
contribution of the parton j to the proton spin. 

The most complete global fit includes all available data taken in spin-dependent deep inelastic 
scattering (DIS), semi-inclusive DIS (SIDIS) with identified pions and kaons, and proton-proton 
collisions. They allow us to extract sets of pPDFs consistently at next-to-leading order (NLO) in 
the strong coupling constant along with estimates of their uncertainties [21 [3]. 

Unlike unpolarized PDF fits, where a clear separation of different quark flavours is possible 
by the use of inclusive charged- current DIS data, differences in polarized quark and anti-quark 
distributions are determined exclusively from SIDIS data and hence require knowledge of the 
hadronization mechanism, encoded in non-perturbative fragmentation functions (FFs). Pion FFs 
are rather well known, but uncertainties for kaon FFs are much larger. Significant progress on the 
quality of fits of FFs is expected once data from B factories and the LHC become available. 

Despite the impressive experimental and theoretical progress made in the last years, many 
fundamental questions related to the proton spin structure still remain unanswered. One of the 
main problems is that present fixed-target experiments suffer from their very limited kinematic 
coverage in x and Q^. The kinematic range could be extended in an Electron-Ion Collider (EIC), 
whose set-up is being currently considered [1], and thus SIDIS measurements would allow us 
to extract Am, Au, Ad, Ad, As, and As with much higher precision. Furthermore, this new 
collider would present the opportunity to perform DIS and SIDIS measurements via charged and 
neutral electroweak currents, which permits to access polarized electroweak structure functions 
[5] that depend on various combinations of polarized quark PDFs and provide an effective way of 
disentangling different quark flavours. 

Experimentally, one has access to the asymmetries, which depend on both the polarized an 
unpolarized structure functions {gi and Fj, respectively). The unpolarized SIDIS structure func- 
tions are well known at next-to-leading order (NLO) in perturbative Quantum-Chromodynamics 
(QCD) [6l[7], as well as the electromagnetic polarized ones [7J. However, polarized SIDIS structure 
functions for charged current and Z boson exchange are only known to leading-order accuracy 
(LO) [H|. In general, a LO calculation only captures the main features but does not provide a 
quantitative description of the process. It is then necessary to know the NLO QCD corrections 
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to both the unpolarized and polarized structure functions in order to extract rehable information 
on the parton distribution functions. 

In this work we present a NLO computation of both polarized and unpolarized SIDIS structure 
functions. We focus on the electroweak ones, which are particularly useful to achieve a full flavour 
separation. The paper is organized as follows. In section [2] we establish the definition of electroweak 
(un)polarized structure functions and present their LO expressions. In section [3] we explain the 
main features of the computation of the structure functions at NLO, focusing on the new results 
for the polarized case. In section H] we present some phenomenological results and analyse the 
relevance of the NLO description of polarized SIDIS. Finally, we present our conclusions in section 

E 



2 SIDIS electroweak structure functions 

In this section we establish the definition of the electroweak structure functions at lowest order 
in perturbation theory. We focus, in this and the following sections, on the semi-inclusive case, 
leaving some comments on the totally inclusive DIS to Appendix [Bl For more details on DIS 
notation, we refer to [9]. 

In lowest-order perturbation theory of electroweak interactions, the cross section for the scat- 
tering of polarized leptons on polarized nucleons and the consequent observation of a hadron H 
in the final state, can be expressed in terms of the product of a leptonic and a hadronic tensor as 

d." ^-V-'J2,,L^^W:^', (2) 



dx dy dz Q 
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where x and y denote the usual DIS variables, — = = Sxy, x = Q'^/{2P-q) (g being 
the electroweak current four-momentum, and S the center-of-mass energy squared of the lepton- 
nucleon system) and z = Ph ■ P/P ■ q the scaling variable representing the momentum fraction 
taken by the hadron H 0. For neutral-current (NC) processes, the sum runs over j = 'j, Z and 
7Z, representing photon and Z exchange and the interference between them respectively, while 
for charged-current (CC) processes the interaction occurs only via the exchange of a. W boson 
(j = W)- The tensor L^^ is associated with the coupling of the exchanged boson to the leptons, 
and the hadronic tensor describes the interaction of the appropriate electroweak currents with 
the target nucleon and the subsequent hadronization of partons into H. The factors r]j denote the 
ratios of the corresponding propagators and couplings to the photon propagator and coupling as 

The unpolarized and polarized SIDIS structure functions {F^^ and gf respectively) are defined 



^We point out that we concentrate here in the current fragmentation region, and thus cuts over the z variable 
must be imposed (typically, z > 0.1). For a description of this process in the target fragmentation region, see Ref. 

m- 
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in terms of the hadronic tensor 



P-q 



2P-q 



+ 



P ■ q 



+ 



P-q 

S-q 
P-q 



g^'\x,z,Q'') 



P-q 



(4) 



where P and S" denote the nucleon momentum and spin four-vectors respectively and P, S are 



P -P -^a S -S -^a 



(5) 



The spin-averaged SIDIS cross section for e^iV scattering and the subsequent production of a 
hadron H in the current fragmentation region, for 3> (where M is the mass of the nucleon), 
is given in terms of the unpolarized structure functions by 



da 



HA 



2710^ 



dx dydz xy 



[1 + (1 - y)'] 2x Pf T [1 - (1 - y)'] xF^^' + (1 - 2 P, 



?H,i 
L 



(6) 



H 



where i corresponds to NC or CC. The longitudinal structure function is defined as P^ 
2xFf^ and vanishes at lowest order according to the Callan-Gross relation [TT]. 77^*^ = 1, while 
T^cc _ (1-1-^)2^^ defined in Eq. ([3])), where A = ±1 represents the electron/positron helicity. 
(For incoming neutrinos, it is rp^ = Arjw instead.) Here and in what follows, the sign ± refers to 
the lepton charge. 



The NC structure functions can be obtained as the sum of the photon, Z, and interference 
contributions: 



^1 



H,NC 



F^£ - [gy ± \g\) V,z F^£^ + {g^;' + g^/±2\g^y g\) F^X 



■,H,Z 



and 



xP 



H,NC 



{g'X ± A g^y) r^,z x Pg^'^^ + [2g^yg^^±X {g^y' + g^')] r^z x P 



?H,Z 
3 ' 



(7) 
(8) 



with g^ = -^ + 2 sin^ 6*^/, g^ 
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For the case of a polarized target, the difference Act of cross sections for the two nucleon 
hehcity states is 



dAcr 



HA 



dx dy dz xy 



HA 



[l + (l-y)^]x gf'' ±[l-il-yy]x g^'' + (1 - y) g^ 



(9) 



where again i corresponds to NC or CC and where g^ = g^ — 2xg^ . Like P^, the latter vanishes 
at leading order. The NC spin dependent structure functions are 
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H,NC 



±A g^'^ - {g\ ± Xgl.) v,z gf''^ + (2 g'y g\ ± Kg^;' + g''/)) Vz ' 



H^Z 



H,Z 



gfr = -{gt^±xgl)v,zg"f + {g'v' + gr±^>^gvgl)vzg', 



H,Z 
L ■ 



(10) 



In the quark-parton model, at leading-order, contributions to the SIDIS structure functions 
Fi and gi can be expressed in terms of the quark (p)PDFs and the hadron FFs. In the NC case, 
these functions are 



-^1 ' -^1 ' -^1 



1 



E K 2 e, gl' + g\^\ (g + g ) 



-^3 ; -^3 ; -'^3 



Y,[^,2e,g\,2g^yg\] {qD^-qDf] 



H,-y H,-iZ H,Z 
9l \ 9l , 9l 



^^^A, ^7^^?^] (AgDf-AgDf). 



(11) 



where is the fractional electric charge of the quark, = ±1 — 2eg sin^ 6'vi/, and (7^ = ±1, with 
the + sign for up-type quarks and the — sign for down-type quarks. In the CC case, since the W 
boson interacts only with certain flavours, we have (assuming four active flavours): 



H,W- 
H,W- 
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uD^ + dD^ + sD^ + cDf, 
2iuDj^-dD^--sD^ + cD^), 
AuD^ + AdD^ + AsD^ + AcD^, 
-Au Dl^ + Ad + As - Ac 



(12) 



For exchange, one should replace u -{^ d and s ^ c. 



3 Next-to-leading order 

Assuming factorization, the polarized SIDIS structure functions (/f^^ at a factorization scale hf 
can be expressed as convolutions of non-perturbative pPDFs Afk{x,fj,p) and FFs Df{z,fip) with 
short-distance coefficients ACC^^^i^x, z, fijp), which can be evaluated in perturbation theory. At 
next-to-leading-order (NLO) in the strong coupling constant as these read 

Qa,qb ^ 

+Df ® ACf ® Ag, + ® ACf ® Ag'^ [x, z, Q^) } . (13) 

where we fix all scales equal to Q. The sum runs over all contributing partonic channels. That 
means, for = 7, 7Z, Z, 

Qa = Qb = u,d,s,c,u,d,s,c (14) 
if we consider four active fiavours, while for V = W~ only 

ga(g6) = u{d) , d{u) , s{c) , c{s) (15) 



standing for 7, Z, 7Z or W^. 
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are considered. For V = W^, one should replace u ■<r^ d and s -H- c. The factor A^^*'" takes the value 
A^^" = 1 for every g^, while it reads A.^^" = A^^" = 1 if Q'a refers to a quark and A^g" = A^f" = — 1 
for Qa representing an anti-quark. In Eq. f|T3|) ® denotes the usual convolution: 

{D®C®f){x,z,Q')= f—D{uj)C{-A.Q')f{y). (16) 

Jx y Jz ^ y ^ 



According to equations ffTTj) and f|T2|) . which represent the leading-order term of Eq. f fT3|) . the 

(17) 



coefficients ACY can be written as 



with 



AC. 



A 



AC. 



AC. 



y 

ol 







Ay 



^ -0 

9 



X 



w 



-A 



w 



In order to calculate the NLO coefficients AC/*^'^, we must take into account the one-loop 
corrections to the partonic process V + qa ^ qt, the real emission V + qa ^ qt + g and the box 
contribution V + g — )• qt + qb- The first two channels contribute to the case with jk = qq, and 
the last one to jTc = qg. The coefficient ACf is obtained from the real emission, considering that 
the gluon is the hadronizing parton. At the intermediate stages of the computation divergences 
appear. In order to regularize them we use dimensional regularization [121 [13], i-e., we work in a 
(i-dimensional space, with d = 4 — 2e. All quarks are considered massless. 

Once the matrix element of each channel jk has been computed, one must obtain the spin 
dependent amplitude 

„ 1 r 9. 9 1 

(19) 

defined in terms of the amplitudes for partons whose polarization is parallel(+) or anti-parallel(— ) 
to that of the target. The calculation of each term in Eq. f lT^ requires projection onto definite 
helicity states of the incoming particles. It is then necessary to make use of the relations 
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u{p, \)u{p, A) 



1 + A75 



for incoming quarks with helicity A (a similar expression is obtained for anti-quarks) and 

PaVfS+PfS Va 



£a{p, Ag) e}{p, Xg 



1 

2fW) 



-9al3 + 



p.T] 



1 pPf]" 

2 p.T] 



(20) 



(21) 



for incoming gluons with helicity A^, where t] is an arbitrary light-like momentum, provided that 
p.T] 7^ 0. The terms independent of A and \g in equations ( 120|) and (I2T]) respectively contribute 
only to the unpolarized amplitude (since they cancel out when subtracting the two terms in Eq. 
( !T9|) ). In the last case, the averaging of gluon spins in d dimensions should be performed by 
dividing by the d — 2 = 2(1 — e) possible spin orientations, as has been made explicit in Eq. fl^ . 
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Once the NLO A M^^ has been computed, we can obtain each one of the coefficients AC/ 



ky 



as the finite part of the partonic structure function, defined by 



1 

47r 



dT Pr A M^' 



(22) 



where dT is the dimensional phase-space and the projectors P^^" are 



PI 



-I e 



2p ■ g' 



Ax" 



^5 



2(1-6) 



(23) 



The functions gj^ contain coUinear divergences, that appear as poles in e (at NLO, simple poles 
in e). We factorize these divergences using the MS scheme, i.e., removing the quantities 



^f(x,z) = AQ 
r(x,^) = AQ 



h 7s - log(47r) 

e 

-- + 1E- log(47r) 



-- + 1E- log(47r) 



[Af,,{x)Sil-z) + P,,{z)6il-x)], 
AP,g{x)5il-z), 



XI 



(24) 



where Pjk and APj^ are the unpolarized and polarized LO Altarelli-Parisi splitting functions [T] 
and 'Je = 0,5772... is the Euler constant. The quantity Afgg is defined below in Eq. fl^ . The 
finite functions obtained after factorization are the coefficients AC/*^. 

The main feature of the calculation described above is the correct use of 75 and the Levi-Civita 
tensor appearing in equations (1201) and fl^Tl) . which is not straightforward in d ^ 4 dimensions. 
For our calculations we use the original prescription of 't Hooft and Veltman [13], afterwards 
systematized by Breitenlohner and Maison [T5] (HVBM scheme), which is the most reliable and 
consistent scheme [151 [161 . 



(25) 



In the HVBM scheme explicit definitions for 75 and e^upa are given. In particular 



75 



4! 



and the e-tensor is regarded as a genuinely four- dimensional object with its components vanishing 
in all unphysical dimensions. The d- dimensional Minkowski space is then explicitly divided into 
two subspaces, a four- dimensional one and a (c? — 4)-dimensional one, each of them equipped with 
its metric tensor, g and g respectively. The Dirac matrices are also split into a four-dimensional 
and a. {d — 4)-dimensional part: 



7/^ T/i ~l~ 7/i 

Each part satisfies the usual anticommutation relation. 



(26) 



As defined above, the 75 anticommutes with 7^ but commutes with 7^. As a result, besides 
the usual scalar products p' -p' {p' being the moment of the outgoing parton), products p' -p' show 
up in calculations, with p' the [d — 4)-dimensional component of the momentum p' . Every time 
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this type of products appear, the two-particle phase space in Eq. ( 122]) must be used in the less 
integrated form [T7] 

dT = ^ {Any -— ^ ^\{p'') ip'Y'^^^- (27) 

" 2 

For those terms which are independent of p' the integration is trivial and one recovers the result 
of [6]. 

When computing the amplitudes, 75 matrices are present not only due to the helicity projectors, 
but also because of the weak interaction vertex. With the definition of 75 given in the HVBM 
scheme, it is important to consistently define the couplings with the quiral fields. It is shown in 
Ref. [IH] that the correct W~ vertex is obtained through the symmetrization 

7^^(1-75) -> ^(l + 75)7^(l-75)=r(l-75) . (28) 
Analogously, the Z vertex is obtained by 

I'ig'v - A 75) = 1^ {g'v - 9a) + 9\ r (1 - 75) -> r {g'v - A) + A T (i - 75) • (29) 



Finally, an important property of the HVBM prescription for 75 is that it leads to helicity 
non-conservation at the qqg vertex in d dimensions, expressed by a non-vanishing difference of 
unpolarized and polarized d-dimensional LO quark-to-quark splitting functions. 



AP 



d=4-2€, 
qq 



X] 



Pi 



=4-2e. 



m 



X] 



4C 



'i? e (1 — x) 



(30) 



As it is discussed in [19], this result entails some disagreeable consequences, such as non-conser- 
vation of the flavour non-singlet axial current and an incorrect result for the 0{as) correction to 
the Bj0rken sum rule. It is then customary to slightly modify the definition of the MS scheme in 
the polarized case and define the quantity 



A/, 



X 



APgg{x)+ACpe{l 



X] 



(31) 



to be removed in Eq. 

We present now the NLO coefficients AC/'^. All the traces have been computed with the 
program TRACER [20j, which masters most of the intricacies of HVBM scheme. All results are 
given in the MS scheme, as defined above. Here, n^p = = fip is made, with fi'^ and fip^ the 
initial state and final state factorization scales respectively. For the sake of brevity we suppress 
in the following results the argument {x, z) of the coefficient functions. For i = 1, 



ACf'^ 
ACf'^ 



y 



\lCw2{l-x){l 



ACf 



gqy 



-AP,,(a;) 
igqy 



V 

— z 



AP,,{x) log ( ^ ^ 

^ fip X 



2(1 -x) 



1 



Cf'" - X\,Cf2{1-x)z; 



(32) 

(33) 
(34) 
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for i = 5 



At/ 



qgy 



AP,,(x)2 



1-z 



V 



Cf'^; (35) 



and finally, the longitudinal coefficients are 



ACf = ^ ACf = 0, ACf = ^ 

Ay A, 



(36) 



The functions are the unpolarized coefficients and are shown in Appendix |X1 The factors A^ 
and Xy are given in Eq. (fT9|) . The quantity AP^^ is defined for simplicity as 



AP,g(x) = 2AP,3(x) = 2x-l. 



Finally, the coefficients ACl^ can be obtained as 



ACf = ACf + 2x ACf 



(37) 



(3^ 



We note that the results for the electromagnetic case [ACl^'^) are in agreement with those of 
Ref. [7]. 



Finally, we make some comments on the unpolarized case. It is convenient to define 

(Jl,J2,3^3) = (2Pl,P2/x,P3) 

and 

and thus, the structure functions can be written at NLO as 



Jf '^(x, z, Q2) = e I ClDl{z, Q') g.(x, Q') + 



2ti 



qq,V 



+Df ® Cf'^ ® + ® (^g]{x,z,Q')] 



with 



cY = A^ 



V 1 '-^03 



(39) 
(40) 



(41) 



(42) 



The factor ^^^^ takes the value ^J'" = 1 for every g^, while it is = ^f^" = 1 if refers to a quark 
and ^1" = = — 1 for g^ representing an anti-quark. 

The computation of these coefficients is similar to the polarized ones. In this case, the spin- 
averaged amplitude must be used, and the projectors are 



>x 



Q 



2 p'p", 



Pi 



1-e 



1 

2' 



P 



I e 



fii/pa 



QpPa 

p-q' 



(43) 



The same care as in the polarized case must be taken when dealing with the 75 matrices present 



in the weak vertices. All NLO MS coefficient functions Cl^'^ are collected in Appendix \^ 
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Figure 1: CC spin dependent SIDIS structure functions ' , 9^ ' , and ' /(2x), at 
= 100 GeV^ for the z-bin z = 0.2 — 0.4. The dashed hues show the LO results (the one for 
gl /(2a;) coincides with that for g'^ ), while the solid curves are NLO. For comparison, we 
also show the electromagnetic g^ 



4 Structure functions at an EIC 



In this section we analyse the relevance of the NLO corrections to the SIDIS structure functions 
we have computed. We focus on the CC case, i.e. the interaction via a W~ boson, since it allow 
us to achieve a full flavour separation. These functions will be accessible in a future EIC, whose 
characteristics can be found in Ref . ^ . We will study the behaviour of the structure functions for 
an energy scale of = loo GeV^ for which the x-range foreseen is approximately 4.10 ^ < x < 1 
(see Fig. 7.16 of Ref. We will consider four bins in the hadronic z variable, with 0.1 < z < 0.8 
and rely on the DSS fragmentation functions set of [2T] . 

In Fig. [U we show the spin dependent structure functions , g^ , and gj /(2x), at 

Q2 = 100 GeV^ for the bin z = 0.2-0.4, using the DSSV set of pPDFs |3]. Results are shown both 
at LO (dashed) and NLO (solid). It is important to mention that NLO pPDFs are used also at 
LO, in order to pick up only the effect of the corrections introduced by the NLO coefficients. One 
observes that the NLO results differ from the LO description, particularly in the small-x range. At 
X = 0.002, for instance, the discrepancies are of about 34% for gi and 15% for g^/ (2x) (these being 
larger than those for g^ for all x). Thus, the NLO description is crucial for a precise extraction of 
the pPDFs in the small-x region, which is particularly interesting since it could be measured at 
the EIC with unprecedented precision. For comparison, we also show the electromagnetic gi , 
for which the NLO corrections are more moderated (below 15% in the mentioned x value). 

In Fig. [2] we present the prediction for the NLO SIDIS structure functions gi'^ with H = 
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Figure 3: CC spin dependent structure function g^'^ for H = tc'^,K'^ at = 100 GeV^ for 
four different z-bins. 

tt'^jK^ for the four z-bins considered. In all cases, the DSSV set of pPDFs where used. We can 
understand the behaviour of these structure functions, at least qualitatively, by considering the 
behaviour of pPDFs and FFs, taking into account that only the four channels of Eq. ( fTSl) must 
be considered, as is discussed next. 
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Figure 4: Left: CC and electromagnetic spin dependent SIDIS structure functions ' (green) 
and gi (red) for z = 0.2 — 0.4. We present the results obtained using DSSV pPDFs set (solid) 
and all pPDFs of DSSV set except for As, taken from DNS set (dashed). Right: The same, for 
CC spin dependent SIDIS structure functions gi for z = 0.2 — 0.4 (green) and DIS structure 
function g^ (violet). 



We focus now on the vr" case. The production of this hadron in the final state is mostly due 
to the hadronization of a quark d (between 30% and 46%) or a quark u (between 28% and 41%). 
This means that the quark present in the initial state must be a quark u or d respectively. Thus, 
the behaviour of the structure function is dominated by that of the Au and Ad pPDFs, the first of 
them being positive for all x-range, and the last one negative. Given that the grater contribution 
is that of the u channel, the vr" structure function gl is always positive. For 7r+, however, at 
low X contributions from Ad and As (negative both) are larger and so g^ becomes negative. 

On the other hand, given the magnitude of the fragmentation functions in DSS se t§, the K- 
production for low z is dominated by the hadronization of a quark c (around 60%), meaning that 
the process is initiated by an s quark. Thus, the shape of g^ is practically that of the As 
distribution, with a sign change in the DSSV set around x ~ 0.02. For 0.4 < z < 0.6 the c and u 
hadronization probabilities are comparable, and for 0.6 < z < 0.8 the last one is even larger, but 
cancellations between Am and Ad occur such that the behaviour of g^ remains very similar 
to that just described. 

In Fig. [3] we show the prediction for the SIDIS structure functions g^'^ at NLO. By looking 
at Eq. (fT2l) we can note that the main difference for these structure functions with respect to the 
previous discussion relies in the fact that these ones depend on (—Am), Ad and As, all of them 
basically negative for all x-range (except for As that is positive for large x, but that effect is not 
noticeable). That explains why gl^''^ is negative for all hadrons considered {H = 7r^,i^^) and 
for all X values. 

§The same being valid for other sets of FFs |22) . 
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Figure 5: Left: CC and electromagnetic spin dependent structure functions g'l ' (green) and 
gl ''^ (red) respectively, calculated using DSSV pPDFs set (solid) and DSSV set except for Am, 
which is taken from DNS set (dashed). In both cases, we use z = 0.1 — 0.2. Right: The same for 
CC spin dependent SIDIS structure functions g'l for 2; = 0.1 — 0.2 (red) and DIS structure 
function g^ (violet). 

As it was explained above, the g^~'^~ structure function behaviour is closely related to the 
As distribution. Thus, one can expect it to be sensitive to a small change in that pPDF. We show 
in the left hand side of Fig. IHthe CC and electromagnetic structure functions g^ (green) and 
g^ (red) calculated with the DSSV set of pPDFs (solid) and the same result after modifying 
only the As distribution according to the one in the DNS set [23] (dashed), as a way to emphasize 
the sensitivity of the observable on that flavour. We note large discrepancies between both sets 
of pPDFs in the CC case, and even observe a region in which they have opposite signs, unlike the 
electromagnetic case, for which the difference between both sets is much smaller. In the last case, 
the presence of an s quark in the initial state implies (at LO) the hadronization of an s quark, 
but its FF into a K~ hadron is almost negligible. The sensitivity of gi~'^~ is also stronger than 
that of the DIS structure function gY for some x-ranges as it can be observed in the right hand 
side of Fig. m 

Similarly, one can expect the g^ structure function to be particularly sensitive to a varia- 
tion of Am. We show in the left hand side of Fig. [5] the CC and electromagnetic spin dependent 
structure functions g^ (green) and g\ (red) respectively, calculated using DSSV pPDFs 
set (solid) and DSSV set with the exception of Am, taken from DNS set (dashed). By comparing 
these two plots, we can note that in fact the CC SIDIS structure function is much more sensitive 
to Am than the electromagnetic one, and thus better flavour separation can be achieved from such 
measurement. 

We also compare in the right hand side of Fig. [5] the CC spin dependent SIDIS structure 
functions g^ ' for z = 0.2 — 0.4 (red) and DIS structure function g^ (violet) obtained using 
DSSV pPDFs set (solid) and DSSV set with the exception of Am, taken from DNS set (dashed). 
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The sensitivity of these functions on Am are similar, but one can still profit from the larger 
structure function in the SIDIS case. 



5 Conclusions 

We have computed the next-to-leading order QCD corrections to the complete set of polarized 
semi-inclusive electroweak structure functions in a consistent scheme. 

We have performed a phenomenological study of possible measurements in an EIC. The NLO 
corrections are found to be important in the small-x range, the kinematical region of particular 
interest at a high energy collider. 

We analysed the sensitivity of SIDIS CC structure functions on different pPDFs, and showed 
that they provide an excellent tool to disentangle the full set of flavour spin dependent distributions 
with an unprecedent precision by including the electroweak polarized SIDIS structure functions. 
In particular, we explicitly show that the structure function g^~'^~ is sensitive to the very poorly 
known As distribution while the function gi~'^~ is particularly sensitive to Au. 
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A Unpolarized SIDIS Coefficient Functions 

Here we list all unpolarized coefficients C/'^'^ for SIDIS as introduced in Section [31 To keep the 
expressions as short as possible we define 

= MO = iiflog«. (44) 

In what follows, we always suppress the argument (x, z) of the coefficient functions. All results 
presented here are given in the MS scheme. 
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The unpolarized coefficients C/'^'^ for SIDIS are 



CI 



qqy 



Xy Cf 



-8 6{l-x)6{l~ z) 



+S{l-x) 

+S{1 - z) 
1 

+2 



P,,{z) log + L,{z) + L^iz) + {l-z) 

P,,{x) log + - L2{x) + (1 - x) 



1 + 2 



l+x 



{1-X)+{1-Z)+ (1-X)+ (1-2) + 



+ 2(1 + X2) 



CI 



qgy 



X 



Cf 



V 



2 

+ Pgg{x 
Xy Cp 



5{l-z) 



Pqg{x) log 



Q^ l-X 

2 



+ 2x{l -x) 



(l-2)_ 



P,,(2) ( 5(1 - x) log ( ^2(1 - z)) + 



+2(5(1 -x) +2(1 + X-X2) 



l + x 



(45) 



(46) 



(47) 



(jqqy 
(jqgy 
(jmy 



Ay 4 Cpxz, 
A|:^4x(l -x), 
A^4C^x(l -2); 



(48) 
(49) 
(50) 



fjqqy 



(jgqy 



CI 



qqy 



C^2(l-x)(l-2) 



V 



2P,,(x) 



1-2 



-A^ 



C^2(l -x)2 



(51) 
(52) 
(53) 



We note that all our results in Eqs. ( H5ll - (J53l) are in agreement with Refs. [HI [7] 



B Unpolarized and polarized DIS Coefficient Functions 



For the sake of completeness, we present the NLO unpolarized and polarized DIS strucutre func- 
tions, py and gY respectively. These, can be expressed as convolutions of non-perturbative 
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(p)PDFs with short-distance coefficients {Cl'^ ^{x, Q^) in the unpolarized case and ACf 5 ^(a;, Q 
in the polarized one). At NLO, these read 



la 



2-K 



Cf ®qa + Cf ®g\ (x,Q2) 



(54) 



la 



2tx 



g,v 



All the relevant factors have been introduced in section [31 



The unpolarized coefficients are 

cf 



^ — XyCp2,X, 



XI. Cf <ilog(^)P,,(x) + 5(l-x) 



9 _ TT^ 

2 3 



3 1 
+Li(x) - L2(x) - -- 



2(1-0;) 



+ 3 



Cf 



Cf 



CI 



q,V 



Cf(I-x) 



Cf 



CI 



V 



V 



XZ2x(l-x) 



Ai/ „ \ Pqg{x 



log 



QH-x 

2 



+ 2x(l-x) L 



0; 



and the polarized ones, 



ACf 



ACf 



ACf 



Ay 



cr-A|;c^(i-x), 



AC 



■g,v 



} 




•gy 



AC{ 



ACf 



X 



AP, 



19 



log 



2 



2(1 -x) 







with Li{x), L2{x) and Pjk defined in Eq. (jUj) and APj^ defined in Eq. (1571) . 
All our results are in agreement with Refs. [HI HTJ EH [25] . 



(55) 



(56) 



(57) 
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